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APPLICATIONS  TO  THE  THEORY  OP  THE  SHOCK  TUBE 

ABSTRACT 

The  theory  of  plane  shook  and  adiabatio  waves  is  presented  in  an 
easily  derived  manner  together  with  sufficient  background  material  to 
enable  the  novioe  in  the  field  to  grasp  the  fundamentals  required  for 
further  study.  The  application  of  the  basic  theory  to  the  shook  tube 
as  a  research  instrument  is  given  together  with  some  experimental  re¬ 
sults  to  illustrate  the  calculations.  Certain  conceptions  of  energy 
and  its  relation  to  the  impulse  in  a  shook  wave  are  presented  in.--.£' 
manner  not  used  in  the  literature.  , ^ 
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INTRODUCTION 


The  theory  of  shook  waves  has  beon  oovered  in  many  exoellent  papers 
by  British  and  Amorioan  authors,  notably  G.  I.  Taylor  and  W.  G.  Penney 
of  England,  and  J»  von  Neumann,  J8  Kirkwood,  G.  KiBtiakowsky  and 
E.  B.  Wilson,  S.  Chandrasekhar,  and  others  of  this  country.  It  is  the 
purpose  of  this  paper  to  present  some  of  their  results  in  an  easily 
derived  manner  together  with  certain  background  material  which  would 
make  it  possible  for  the  uninitiated  quiokly  to  gain  a  oertain  per¬ 
spective  in  the  field  with  the  objeot  of  applying  this  baokground  to 
applications  of  the  shook  tube  as  a  research  instrument.  Certain  oon- 
oeptions  of  energy  and  its  relation  to  the  impulse  in  a  shook  wave  are 
presented  in  a  manner  not  used  in  the  literature,  and  experimental  re¬ 
sults  are  quoted  to  show  the  validity  of  oertain  assumptions  neoessary 
for  such  relationships  to  hold. 

The  use  of  a  bursting  diaphragm  in  a  tube  as  a  method  of  producing 
shook  waves  is  covered  in  ft  phenomenological  report  by  W.  Payman  and 
W.  Co  P.  Shepherd  in  1941.  They  attribute  its  earliest  conception  to 
P.  Viello  in  1899.  A  calculation  of  the  shook  wave  pressures  to  be 
expeoted  thereby  was  oarried  out  by  A,  H.  Taub  in  1942.  Certain  experi¬ 
mental  work  on  the  tube  as  a  primary  standard  of  blast  wave  pressures' 
was  done  by  G.  T.  Reynolds,  and  extended  in  a  series  of  quite  aoourate 
measurements  by  W,  T.  Read  who  found  that  the  experimental  pressures 
varied  about  6  peroent  from  those  calculated  from  tube  theory. 

L«  G.  Smith  has  used  the  tube  as  an  aid  in  studying  experimentally  the 
reflection  of  shock  waves  at  oblique  incidence.  It  has  been  used  as 
an  instrument  in  the  study  of  bursting  diaphragms  from  incident  shook 
waves  and  as  an  adjunot  to  the  development  of  piezo-electrio  pressure 
gauges.  Certain  proposals  have  been  made  for  using  the  tube  on  a 
fairly  large  scale  for  tests  on  land  mines  and  other  devices  exposed  to 
explosivo  blast  in  an  effort  to  reproduce  field  results  in  the  labora¬ 
tory. 


It  is  the  earnest  hope  that  this  paper  may  provide  sufficient  back¬ 
ground  so  that  persons  who  have  not  previously  used  the  tube  may  readily 
understand  in  a  quantitative  way  the  phenomena  that  may  be  observed, 

I.  THE  PROPAGATION  OF  FINITE  ADIABATIC  WAVES  IN  A  TUBE 

The  derivation  of  the  properties  of  one -dimensional  finite  waves 
oan  be  oarried  out  quite  readily  by  oonsidering  a  slice  of  gas  in  a  tube 
bounded  by  planes  at  x  and  at  x  +  dx  as  shown  below. 


u 


$P 


x  +  dx 


The  gas  in  the  thin  slioe  is  assumed  to  have  a  density  and  a 
particle  velooity  u,  both  of  whioh  are  functions  of  time.  Then  one  oan 
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whore  u  «=  f  (t) 


write  for  tho  foroes  on  the  alioe 

/*&  *  ^(u)  =  -  $  P  where  u  «=  f(t) 

du  «  Su  a.  A*  flU  £u  .  $P  dP 


<s»x 


fix  sx 


So  from  equation  (1*1)  we  have  the  familiar  equation  of 
a  gaa  in  a  tube. 

3u  .  „  au  .  1  <3P  rt 
c?t  +  a*  ax  ® 

The  equation  of  continuity  ia 

-  3  /  „ ^ 

ax' 


but  ainoe  jO  xz  f  (u) 

and  4£-=4£.3" 
ax  du  ax 


W  -  -ST 
f£+"-§£v<t^o 

we  have,  after  making  indioated  changes  of  variable, 

tHI- 

Cancelling  the  equation  of  continuity  may  be  written  as 
3  u  ...3u  .  ~  du  su  _  au  .  ,  . .  _  du  %  3u  ~ 

stt*  f  JT-si'ir*  <“+/>3^->-s*’=° 

Since  P  =  f(/&)  and  /»  =  f(u)  we  have  that 

d 


3F  dP  a/O 
«?x  HJo  3  x 


and  a/» 
3  x 


au 

a  x 


so  that 


J>jP  dP  d£  j£u 
3x  ”  d/9  au  3x 


Substituting  this  expression  into  the  equation  of  motion  (l.J 
that  it  may  be  written  as 

3u  .  „  3  u  1  dP  d/O  3u  _  3u  ,  .  1  dP  &/o  3u 

These  two  equations  (l.2a  and  1, 3a)  may  be  made  consistent  il 

-  du  1  dP  d fi  /du  /  _  1  dP 

«*  %o>  ~~pZT? 


(I.D 

motion  for 

(1.2) 

(1.3) 

(I.S*) 

)i  we  see 
0  (I. 2a) 

(1.0 
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then 


and  therefore 


d/»  -/° 


1  /dF 


/© 


(1.40 

(1.5) 


This  is  the  expression  for 'the  particle  velocity  in  terms  of  the 
pressure  and  density  for  waves  which  travel  in  one  direction. 

Returning  to  the  equation  (l«2a)  of  motion  we  see  that 

from  equation  (l»4a) 

so  that  we  have 

■§!£  +  (u  +  =  0  whioh  equals 

■fr 4  (u  t  fxp*  §£  *  0 


(1.6) 


For  an  adiabatio  compression  or  expansion  — -? 

/°r 


P  _ 


dP _ r-i  _  p 


rP 

/• 


yp  p 

but  a  ’  where  a  is  the  looal  velocity  of  sound  in  the  medium  then 


/dP  _ 


and  the  equation  of  motion  reduces  to 

3u  ,  /  3  u  * 

sT  +  (u  J,  a)  _  =  0 


(I. 6a) 


If  we  wish  to  find  the  velocity  of  a  section  of  the  wave  of  con¬ 
stant  particle  velooity  and  consequently  oonstant  pressure  we  oan  do  so 
by  letting  u  be  a  constant  so  that  du  =••  0 


but  slnoe  u  =  f (x  t) 

9u 

.  dx  _  3 1 

so  that  TC  *  -  "U 

ax 

But  from  (X.6a)  we  see  that 

■5; 


*>  =  -§Tdt+§i,s*  =0 

velooity  of  propagation  of  that 
seotion  of  the  wave 
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so  that  the  velocity  of  propagation  of  a  section  of  the  wave  of  oonstant 
partiole  velocity  is  equal  to  (a  +  u)« 


v  e  a  +  u 


(1.7) 


The  sign  of  u  is  positive  if  the  particle  velocity  has  the  same 
direction  as  the  wave  propagated  in  the  medium  and  is  negative  if  the 
two  are  in  opposite  directions.  The  looal  velocity  of  sound  a  will  he 
a  function  of  the  pressure  P  as  will' the  particle  velooity  u. 

If  the  wave  is  an  adiabatio  compression  (not  a  shook  wave)  then  we 
may  evaluate  the  partiole  velooity  u  and  the  looal  velooity  of  sound  as 

followst  Assume  that  P  «=  P  (^)  for  an  adiabatio  pressure  ohaage 

°  /o 


4%  1  °  V 


a?  ryw 

,  rr*  r, 

ro 

\ 

r  rP* 

\}~t 


.  2  /  fo. 


rP  r-j 


P  * 

f  o 


sr{/*  -/“« 


f—  y-1  -1 

'  ‘] 


(1.8) 


r  p0  ^ 

but  (—75—)  =  a  the  velooity  of  sound  in  the  mediun  into  whioh  the 

'O  0 

wave  advanoes.  T'-l 


So  ’that 

- 

n 0  0 

oonsequently 


--lifts 


—  ""A  — I 

(-^")  -  1  (1.1 

L^o 


u  *i^r=3 


4? 


(l.8b) 


partiole  velocity  for  a  oompressional  wave. 


In  the  same  manner  we  find  that  the  partiole  velooity  for  a  rare¬ 
faction  wave  is  7  -1  ~ 

2  p  “sr 

u  -  +  yzi  a0  1  *  (|r*)  (1.8) 


I7TT 


(1.8) 


Iffhere  a  is  again  the  velooity  of  sound  in  the  medium  into  whioh  the 
o 

wave  advanoes. 
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The  local  -velocity  of  sound  a  is  related  to  the  -velocity  of  sound  in 

r  -I 

/°„  r?„  -y  rro  „  '  r 

(jr) 


a  medium  of  pressure  PQ  and  density  /30  as  follows) 

a2  *22«1!2  P  «  'll*  ,  <*)(*) 
^  0  Ai  0  v'v 


so  that 


a  as  a 


* 


2* 

>  < 


dao) 


The  velocity  of  a  section  of  oompressional  wave  of  constant  partiole 
velocity  advancing  into  an  undisturbed  medium  is 


V  * 
uncons 


$ant- 


a  *a  a 


r-i 

(p^+  2 
>1f~'  +  F=r  \ 

o 


If  ^  e  1,4  then 

1/  — ,  a 

"u=oonatant  0 


(I, 10a) 


This  ia  the  velocity  also  for  a  rarefaction  wave  advancing  into  a 
medium  of  pressure  Pq  and  density  /a  Q  and  partiole  velooity  uQ  «  0» 


Prom  equation  (l*10a)  it  can  be  seen  that  the  higher  pressure  parts 
of  a  wavs  will  travel  faster  than  the  lower  pressure  parts  so  that  a 
finite  compreBBional  wave  would  assume  the  shapes  shown  below  in  time 
sequence* 


a'  b*  o1  d* 

A  oompressional  wave  tends  to  assume  a  steeper  slope  at  the  front 
of  the  wave  while  ararefeotion  wave  tends  to  become  less  steep  and  to 
change  shape  in  the  opposite  direction*  The  case  (d)  of  the  compress- 
v ional  wave  can  not  happen  for  the  front  has  a  double  value  of  pressure 
but  ^he  front  doeo  become  vertical  with  the  formation  of  a  shook  wave 
which  dees  not  obey  the  same  equations  as  does  the  adiabatic  waves* 
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The  rarefaction  waves  however  obey  the  adiabatio  equations  at  least  to 
the  trough  of  the  wave*  It  can  be  seen  that  there  Is  a  tendency  to  form 
a  secondary  shook  in  the  tail  of  the  rarefaction  wave*  This  phenomenon 
will  not  be  treated  here.  The  shook  wave  when  formed  obeys  a  set  of 
conditions  embodied  in  the  Rankins -Hug oni ot  equation  whioh  will  be  de¬ 
rived  in  a  suooeeding  section. 

IT.  THE  REFLECTION  AND  TRANSMISSION  OF  ADIABATIC  WAVES  AT  A  DISCONTINUITY 
(a)  Compressions!  Waves. 

The  sequence  of  events  where  an  adiabatio  wave  meets  a  discontinuity 
oan  be  shown  by  the  following  sketoh. 


-v 

ul— 


f 


denser  medium 


PjjJ^-iees 

med 

- __ 


less  dense 
medium 


o«. 


The  velocity  of  sound  on  the  right  of  the  boundary  is  a  end  on  the 
left  of  the  boundary  is  a^«  The  pressure  initially  is  the  same  on  both 

sides  and  is  P^.  .An  adiabatio  wave  moving  to  the  right  with  a  pressure 

Pp  partiole  velooity  Up  and  velocity  of  sound  behind  front  of 

approaches  the  boundary.  Depending  on  the  oonditions  a  wave  of  greater 
or  lesser  amplitude  will  be  transmitted,  while  a  compressional  or  rare- 
faotian  wave  will  be  refleoted  to  the  left.  The  pressure  and  the 
partiole  velocity  behind  the  transmitted  wave  are  Pg  and  Ug,  respectively, 

while  the  pressure  and  the  partiole  velooity  of  the  refleoted  wave  are 
P3  and  u3. 

The  two  oonditions  whioh  must  be  satisfied  at  the  boundary  are 

Pg  =  P3  pressures  equal  (ll.l) 

and  ul  “  u3  *  ug  partiole  velocities  equal  (II.2) 


The  values  of  the  particle  velocities  are  (if*/  **  1.4) 

r»i 1/7  i 
ui c  6ai  Mp)  *  A 


(n.s) 

(n.4) 
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I 


i 


px  V* 

but  »2'c»i  (^-) 

o , 

rP,  *a '  '  px 

•  O  Uj  e  S»x  jjy-)  -  !f-)  J 


(ZI«6) 


(11.5a) 


Than  from  equation  (ll.2)  we  havw 

», X/1 


t  lA 

5mi  I  (yr)  -  1  -  (p~) 
o  o 


*l[< 


PllAl  fP2lA  1 

4  (y~)  J  *=  69  j  (j j-)  *•  Ij  (H*6) 


and  since  Pg  “  P3  from  equation  (iZ.l)  we  have 

1/7 


P2vlA  2  ai  ,P1 


iTTT 

1  o 


a  -  a. 


»  +  a. 


or  in  more  useful  form 

2  a. 


i A  - 


c  +  ax  ‘1 


1/7  +  *  ~  A  p  1/7 


a  +  a2  0 


(11.7) 


(U.8). 


If  the  boundary  is  a  rigid  wall  so  that  Ug  *=  0  then  we  have 

«1  =u3 


?,lA  =  a^A  -  po1A 


(H.») 

(11.10) 

(11.11) 


(b)  Rarefaction  Waves. 


We  show  a  similar  sketch  for  the  sequenoe  of  events  with  a 
slightly  different  notation  to  adapt  it  to  the  notation  in  the  other 
sections. 
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Where  P  and  c  are  the  pressure  and  sound  velooity  on  the  right 
of  the  boundary  ana  F  and  a^  are  the  pressures  and  sound  veloolties 

on  the  left  of  the  boundary.  A  rarefaction  nave  of  pressure  Pr  ^ 

particle  velooity  Up  and  sound  velooity  behind  the  front  of  a,,, 

approaches  from  the  left.  Depending  on  the  values  of  a^  and  o  a  rare* 

faction  or  a  small  compression  nave  is  refleoted*  while  a  rare faot ion 
of  greater  or  lesser  amplitude  is  transmitted.  The  pressure  and 
partlole  velooity  behind  the  transmitted  wave  is  IV  ^  and  Ug  while 

those  behind  the  ref looted  wave  are  Pr  and  Uj  respectively. 


The  conditions  to  be  satisfied  are 


and 


where 


"5 

1/7 


but 


so 


(V) 


*1  [< 


Pr  lA 

<-A 


Prs  *A 


] 


Equating  the  partiole  velocities  we  have  sinoe  Pr 


(Pr)1/7 

2 


2>i 
C  +  a. 


Pr 


1 A  +  till  plA 


(il.la) 

(ll.2a) 

(11.12) 

(11.15) 

(11.14) 


(11.14a) 


c  +  a, 


(11.15) 


Similarly  for  reflection  against  a  rigid  wall  where  ug  e  0  we 


hav©  the  conditions 


Pr  v  2Pr  ^  (11.16) 

3  1 

These  arc  seen  to  bo  identical  in  form  with  those  derived  for  com¬ 
pressions!  waves.  <’> 

ID. 'DERIVATION  OF  THE  RANKINE-HUGONIOT  EQUATIONS  FCR  A  SHOCK  WAVE 

If  a  shock  wave  travels  with  a  velocity  U  into  undisturbed  air 
(air  at  velocity  uQ  *=  0  in  which  the  pressure  is  Po  and  the  density 

is  /O  )  and  if  the  air  behind  the  shook  is  at  a  pressure  P,  density/^, 

and  is  moving  with  a  velocity  u,  then,  by  using  the  fact  that  for  a 

unit  mass  of  air  orossing  the  shock  front  we  must  have  mass,  momentum 

and  energy  conservation,  we  obtain  equations  involving  U,  P,  p  #  v, 

P  ,  p  and  u  whioh  are  known  as  the  Rankine-Hugoniot  equations.  We 
o/o  o 

derive  them  as  follows i 


1 - 

I~JU 

lp . /°  u . 

Po  /°0  u0 

Let  us  oonsider  an  observer  that  moves  with  the  shook  front.  In 
one  second  the  amount  of  matter  that  crossed  (from  the  right)  a  unit 
cross  section  of  the  wave  front  is  Q  U«  This  must  equal  the  amount 

that  gets  away  frctn  the  left  face  of  tue  oross  section  in  the  same 
interval  of  time,  viu.,^0  (U  -  u).  Henoe  we  obtain  the  conservation 
of  mass  equation. 


jO  (U  -  u)  »  /O  U  =  m  conservation  of  mass  (Ill.l) 

The  momentum  of  the  mass  /O  U  is  /o  U2  and  the  momentum  of  the 

2  0  '  0 

j  iss  /0(U  -  u)  is  /0{ U  -  u)  «  The  change  of  momentum  aoross  the  shock 
iront  must  equal  the  force  acting.  This  is  the  difference  of  pressure 
the  two  sides  of  the  front  times  the  cross  section  which  we  have 
taken  as  unity.  Hence  we  obtain  the  conservation  of  momentum  equation 

p  -  p0  -/yj2  ~/>(U  «  u)2 

or 

P  +  />(U-u)2  =P  +/0  u2  observation  of  (II1.2) 

/  x  *  o/o  momentum. 

To  obtain  the  energy  equation  we  need  to  know  the  internal  energy 
of  the  unit  mass  of  gas  (when  it  is  at  a  pressure  PQ  and  density  /©0)» 

This  is  the  work  done  against  external  pressure  when  the  gas  i*s  expanded 
adiabatioally  to  sero  density.  We  call  these  internal  energies  E  and  Eq 
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,  The  work  done  by  pressure  per  unit  urea  per  seoond  on  a  oolvaon  of 
gas  of  unit  oross  seotion  (the  oolusm  extends  through  the  shook  front) 


is 


p  •  u  -  P  (U  -  u) 


This  oust  equal  the  ohange  ih  kinetic  energy  plus  the  ohango  in 
internal  energy  of  the  gas*  The  former  is 


■g-  m  jjTT  -  u)2  -  TJ2^ 


and  the  latter  is 


»(E  -  B0) 


where  E  and  E  are  the  internal  energies  of  the  gas  on  the  left  and 
o 

right  sides  of  the  shook  front  respectively  and  m  is  the  mass  of  air 
crossing  the  unit  oross  seotion  of  the  front  per  second  and  is  given  by 


m.*=  f>Q  U  «=^o(U  -  u) 


Henoe  we  have 


Po  U  -  P(U  -  u)  *=  -j  jo  -  u)2  U2 J  +  »  (E  -  E0) 


dividing  by  m  we  have 
P 


7- 


or 


S-  *  i  «*  *  is.  -  %r  (»  -  »>2  *  *  •  T.lZnUon  of  d»:3) 


energy 


For  an  ideal  gas  the  internal  energy  may  be  calculated  as  follows  1 


P  «=  k^o ' 


where  /O  « i- 


Henoe  E  » P  d(— ■)  =  -  ^  &  “ky^  ^  d^o 


Tr 1 


therefore 


E 


k  = 

1  P 
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end 


E 


o 


1 


Substituting  these  values  for  E  and  Eq  into  equation  (111*3)  we  obtain 


+. 


1  ro 

^  K 


or 

^  +  s  ^  ^  jsr  +  *  [°  *  “] 2 

Now  we  nay  use  these  fundamental  relations  to  derive  the  relation 
between  the  velooity  of  the  shook  wave  and  its  pressure  together  with 
oertain  valuable  information  ooncerning  the  gases  behind  the  shook* 
from  equations  (lll*l)  and  (III*2)  it  follows  that 

P  -  PQ  =  m  [u  -  (U  -  u)J 


henoe 

-  \  (p  -  PQ)  [U  +  (U  -  u)J  «=  |  O  **  u)2  -  U2J 


and  from  the  third  equation  (III. 3)  we  have 

PQU  -  P(U  -  u)  -  J  f(U  -  u)2  -  U2  J  m  m  (E  -  EQ) 

substituting  from  above  we  have 

PQU  -  P(U  -  u)  +  j  (P  -  PQ)  [u  +  (t?  -ji)  ^  «  m  (E  -  Eq) 

which  when  cleared  beoomes 

(P  +  Pq)  U  "  OLrJl).  =  e  -  Eo  but  a  =  qU  *=  />  (U  -  u) 

whioh  then  is 


This  equation  may  be  interpreted  as  saying  that  the  increase  of 
internal  energy  across  the  shook  front  is  due  to  the  work  done  by  the 
mean  pressure  in  performing  the  compression. 

A O 

from  equation  (I II. 4)  we  may  solve  for  ^  whioh  we  do  by 
multiplying  the  expression  by  and  dividing  by  then 
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•jrU  +  {-)  irC-  -  1)  *  (J-  -  *£-) 

o  '  O  O  '  A 


Let  —  t=  y 

and  ib  z 

.0  ' ' 

Than  it  may  ba  written 

|(y  +  l)(x  -  1)  *  —y  (y  -  x) 
the  solution  of  which  is 

x  g^l  *  [r-ily  "  4"  r*ti0  of  d#nBiti«»  (iii,6) 

'  '  o 

From  aquation  (ill.l)  we  have 

/° 

U  -  u  =  X7 

Henoe  aquation  (ill .2)  nay  ba  writted  as 

p  +  ~jsr  &  B  p0  +  /®0  ^  w  o  ^  *  p  -  p0 

u2  = 


or 


P  -  P 


£ -  K  ®  (2-1^) 

4%,  ™ 


^ c(1  "  7s" 

rp0  2 

Let  — -r—  *=  a  where  a  is  the  velooity  of  sound  in  the  undisturbed 
/°0  0  0 
medium. 

^  TJ2  11. 

— g-  s=  (2  ■"■  y)  end  substituting  from  above  we  haws 

*c  '  1  •*  —  ... 

U2 


“ 
a0 


Therefore  the  first  of  the  important  derived  relations  is  the  velocity 
pressure  relationship  for  the  shook  wave 


t  £p-i  +  (^+l)yj  velooity  of  shook  ware  (ill. 6) 
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FJrora  aquation  (lll,l)  we  have 


henoe 


or  n  «s 


u  B  (1  -  -0)  u 

/>o 


v  r 

velocity. 


ft  «  i  -  4°  c  j  -  llZr-ifc 

U  r-1  4  (;T4l)y 

-T24y,('/ ’ll )y  of  Petiole  velooity  to  shook  (III, 7) 


If  o  s, 


rrp 


where  o  la  the  velocity  of  sound  behind  the  shook  front 
2  /°o  P 


o-  _  ^  P 

3  r7> 


o 


2 

o 

3 


y  Lr  -l  V  (r+i)y 


ratio  of  velooities  of  sound  behind  and  in 
front  of  shook  wave.  (I II . 8) 


rVam  equation  (ill. 6)  and  III. 7)  we  can  derive  the  ratio  of  the  partiole 
velocity  behind  the  shook  front  to  the  speed  of  sound  in  front  of  shook 


2  g 

J  grLr---^-(Wl)yj  <««•»> 

These  equations  although  derived  for  a  shook  wave  moving  into  still 
air»  uq  =  0,  will  hold  for  a  shook  wave  moving  into  air  traveling  with 

a  uniform  volooity  uq#  if  we  understand  U  to  mean  the  velocity  of  the 

front  relative  to  the  moving  air. 


If  the  medium  into  whioh  the  shock  waves  travel  is  air  which  has 
a  T  equal  to  1.4  these  equations  may  be  simplified  and  rewritten  using 

Ps 

I  c  -2.  *=  y  -  1 


/°  _  1  4  6  y  _  7  4  6  2 

/»0  6  4  y  i  4  z 

(111,6ft) 

U2  1  4  6  y  ,  .  6 

4  =1  +  "  ' 

(ill, 6a) 

v=  -  rUs 

(III. 7a) 

■j  =  ^rriy)  -  <* +  h  <77+Y,> 

(III.8a) 
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(HI, 9*) 


u2  _  25{y  -  l)2  _  25  s2 

.2  7(1  +  6yS  y(7  T 

*o 

Three  of  these  quantities  aro  plotted  in  the  following  pages 
(Graphs  1,  2,  and  3)  as  funotions  of  the  exoess  pressure  behind  the 
shook  In  pounds  per  square  inoh*  The  exoess  pressure  is  (P  -  PQ) 

where  Pq  is  assumed  to  be  14,7  pounds  for  computational  purposes. 

\  » 
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Velocity  of  Shock  Wave  as  a  Function  of  Excess  Pressure  in  Shock  Wave 


Graph  3 


1- 


17.  HEFLECTION  OF  SHOCK  WAVES 

•  « 

Whan  a  shook  wave  strikes  perpendicularly  on  a  flat  surfaoe,  a 
,  shoojc  wave  is  refleoted.  The  pressure  at  the  surfaoe  jumps  instantane¬ 
ously  ftrom  the  atmospherio  pressure  PQ  to  the  pressure  P  behind  the 
refleoted  shook  wave.  At  a  point  close  in  front  of  the  refleoting 
plane  the  pressure  first  changes  to  F  the  pressure  behind  the  inoident 
shook  wave  and  then  to  P  as  the  refleoted  wave  reaches  it.  The 
relationship  between  y  *=  P  /fc.and  y  e  P/P  ran  be  derived  from  the 
shook  wave  equations.  '* 


Before  refleotion 


^air  velooity 
sound  velooity 
pressure 


After  refleotion 


' 


air  velooity 
sound  velooity 


pressure 

TT  * 
r 


u-* 

0 

P 


Diagram 
showing 
positions 
at  which 
symbols 

*ppiy. 


If  u  is  the  velooity  of  the  air  behind  the  inoident  shook  wave. 
U  the  velooity  of  propagation  in  still  air,  a  the  velooity  of  sound 
in  the  undistrubed  air,  U  the  velooity  of  propagation  of  the  re¬ 
fleoted  wave  and  o  the  velocity  of  sound  in  air  behind  the  inoident 
wave,  the  shook  wave  equations  for  the  inoident  wave  arei 


si  =  ^  [r-l  +  (r+1)jr] 

O 

o2  -  (r  -I)y1 

a  2  ~  y  Lr-i  +  tr+ijyJ 

o 

«L  -  _  2(y  ~ 1)2 

ft  2  "  V£r-i“+  &  +i)yj 


(iv.i) 

(IV. 2) 


(IV.3) 


(xv.4) 


The  refleoted  shook  wave  advanoes  with  a  velooity  u  4  U  relative 
to  air  in  whioh  the  velooity  of  sound  is  o.  The  velooity  of^the  air 
behind  the  refleoted  wave  relative  to  that  in  front  is  u  as  in  the 
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^  - 


inoident  wavej  thus  the  equations  for  the  ref looted  wave  arei 


then 


u 


2(yr  -  l) 


U  +  Ur  r-1  +  ~(7+l)yr 


u  _  2<yr  ■ 

^  -  r+i  +  V-ijyr 


r  r+1  +  (r-l)yr1 
V  |_r-r+"“(r'+T)3rJ 


—?  55  y 

o 


and  in  reflected  wave 

u2  _  2  ^ 

2  =  7  >  -1)  4  (r+l)y, 


(IV.5) 


(IV. 5a) 


(IV. 6) 


(IV.7) 


r 


Here  o  is  the  velocity  of  sound  in  the  air  behind  the  reflected 
wave.  The  relationship  bet/een  y  and  y  is  found  by  eliminating  u, 
a  ,  and  c  between  equations  (IV.3),  (IV^4),  and, (IV.7).  ft-cm  (TV.3), 
and  (IV.4)  we  have 


u 


_  2  (y  -  l  r 

7  '7* 


(IV.  7a) 


Equating  this  to  (IV, 7)  we  have 

(y  -  l) 


(yr  -  i)' 


iy  -  i)  wr 

yjV+1  +  (y-l)yj  "  r-H  (r+l)yr 


(iv.  8) 


If  one  assumes  ^  for  air  equal  to  1,4  then  equation  (IV*8)  becomes 

v2 


.CZ..--.U  -  s 

y(6  +  y} 

This  oan  be  expanded  into 


(yr  -  l) 

1  +  6  yr 


(IV.  8a) 


y2yr2  -  8  yry2  +  6  y  yr2  +  8  y  =  6  yr  +  1 


2..  2 


and  faotored 

8  y(l  -  y  yr)  -  6  yr(l  r  y  yr)  e  1  -  y*Vr 

whioh  equals 


8  y  -  6  yr  =  l  +  yyf 
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(IV.  8b) 


Pressures  for  intermediate  strength  shocks  are  plotted  an  the  next 
page.  (Graph  4) 

V.  THEORY  OF  SHOCK  WAVE  FORMATION  IN  A  TUBE  FROM  A  BURSTING  DIAPHRAGM 

Consider  a  tube  of  constant  cross  section  closed  at  one  end,  with 
a  gas-tight  diaphragm  fixed  at  some  point  in  the  tube  in  suoh  a  way 
that  a  section  of  the  tube  bounded  by  the  closed  end  and  the  diaphragm 
may  be  pumped  up  to  a  pressure  Pc,  while  the  remainder  of  the  tube  re¬ 
mains  at  a  pressure  PQ.  Then  if  tho  diaphragm  is  suddenly  broken  by 
air  pressure  or  other  means  a  shock  v/ave  will  be  formed  in  the  low 
pressure  section  of  the  tube  advancing  along  the  tube  away  from  the 
diaphragm.  At  the  same  time  a  rarefaction  wave  will  be  formed  in  the 
high  pressure  section  of  the  tube  which  will  progress  baok  into  the 
high  pressure  gas  until  it  is  reflected  at  the  closed  end  of  the  com¬ 
pression  chamber.  Meanwhile  the  shock  wave  will  progress  down  the 
tube  until  it  is  reflected  with  either  positive  or  negative  phase  at 
the  other  end  of  the  tube.  The  sequence  of  events  may  be  shown  in  the 
following  way i 
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Pressures  in 


Graph  4 


icidont  and  Reflected  Shock  Waves 
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^.diaphragm 


before  breaking 


at  breaking 


after  breaking 

region  in  which  pressure  is  P,  ve¬ 
locity  of  gas  is  u,  velooity  of  sound 
is  o  and  velooity  of  shook  wave  is  U. 

Rrom  the  Ranklne-Hugonoit  equations  we  have  the  relationship  be¬ 
tween  the  particle  velooity  behind  tho  shook  front  and  the  pressure 
ratio  P/l»o  at  the  shook  front.  This  isi 


u 


(v.l) 


a 


-  P 


y .  r -i  + 

0  /P o 

where  y  =  4—  ,  a  =  If  — —  f  and  u  =  particle  velooity  behind  shook  front 
o  °  *  o 

The  region  behind  the  diaphragm  is  propagating  a  rarefaotion  wave 
beoause  of  the  relief  of  pressure  by  the  bursting  of  the  diaphragm.  The 
velooity  of  the  particles  in  the  rarefaotion  wave  where  the  pressure  is 
P  is  given  by 

_ 

d/0  which  has  been  evaluated  in 

/&  1  equation  (1.9),  and  found  to  be 

Y- 1 

2  _  I  ,  /P  \  Tr 


U  =  + 


=  i  r~=T 


r  ~~ 

i  -  <f. )  ^ 

_  o  _ 


(V.2) 


where  P  i3  the  pressure  of  the  region  into  which  the  wave  advances  and 
a^  is  tSe  velooity  of  sound  in  this  region. 


Now  P  =  P 


P 

FT 


=  V 


so  equation  (V.2)  becomes 


u  =  + 


-  r^T 


r- 1 

~TF  1 


i  -  (y 

o 


(V. 2a) 


These  two  particle  velocities  must  be  equal  (at  the  diaphragm  after 
breaking);  otherwise  a  local  region  of  vacuum  or  high  pressure  will  de¬ 
velop  in  time.  So,  equating  the  two  particle  velocities  we  have 


2 


2ao2  (y  -  1) 

rrr-i  +  ( r^nTyj 


4a, 


( r  -D‘ 
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r-i  2 
Tf 


i  -  (y  p2) 


(V.3) 


Assume  a  =  a-^  which  moans  that  sufficient  time  must  elapse  after 
pumping  up  the  compression  chamber  for  the  gas  to  be  reduood  to  ambient 
temperature*  Extracting  tho  square  root  on  both  sides  we  have  then 

:  .  -i*i 

_/2.(yi?  e  2 

V -  1  +  (y  +  I)yj  *  T 

Assume  y  =1.4,  then 

JL  z  2 


(y  5™) 

O 


*=  X 


?olA 

(y  jr) 


of 


7f 


1  - 


y  -  i 

/  7(1  +  6  y) 


V  7(1  +  6  y) 

Implicit  solution  for  y 


(V.4) 


This  equation  then  gives  the  pressure  in  the  compression  ohamber 
before  breaking  the  diaphragm  necessary  to  establish  a  shook  wave  of 
pressure  P  travelling  down  the  tube.  This  relationship  is  plotted  on 
the  following  page.  (Graph  5).  On  the  second  page  following  is 
plotted  also  the  values  of  exoess  ohamber  pressure  over  atmospheric 
and  the  exoess  shock  wave  pressure  over  atmospheric  in  pounds  per 
square  inch.  (Graph  6). 


The  compression  ohamber  pressure  may  also  be  evaluated  in  terms 
of  a  given  Mach  number  M  behind  the  shock  front. 

Since  ^  =  M  t=  ......  -  . —  -  **■  from  equation  (ill. 9a)  (V.6) 

/7(7  +  6i)  y7(l  +  6  y) 

Equation  (V.4)  may  be  rewritten  as 


y  __  z  +  1 

(1  -  I)'  (1  -  f)? 


(V.6) 


Solving  (V.5)  for  z  in  terms  of  M, 

z  =  0.84  M2  +  1.4  M  /l  +  0.36  M2 


Pc 

_  1  +  0.84  M2  + 

1.4  M  /TT 

0.36  M2 

F* 

0 

(1  * 

uJ 

’V 

M 

P<A 

P  /P 
s'  o 

1.4 

52 

4.25 

1.2 

29.3 

3,3 

1 

16.5 

2.5 

29 

_  „ _ ■ —  •  ’  . 

(V.7) 


In  oaoe  gases  of  different  density  and  y  are  used  in  compression 
chamber  and  expansion  chamber  the  relationship  between  compression 
ohamber pressure  and  shock  strength  will  be  different.^  the  gas  in  the  com¬ 
pression  uiininutir  nas  a  gSufisa  y ^  and  scun.u  velocity  and  the  gas  in 

the  expansion  ohamber  has  a  ganma  and  sound  velocity  a?  then  it  can 
be  shown  that  %  • 


r0 


1  «* 


(y-i)  (7.-1)  T57! 


(v.e) 


For  example,  if  helium  is  used  in  the  compression  ohamber  and  air  in  the 
expansion  ohamber, 

ri  =  1.66 

72  =  l.*o 

ai 

~  «=  1.93 
a2 

At  y  *=  3, 


P 

o 

r 


"f  ■-  1  . - . "aw  ■  ■KT  *6*7 

_  _ .518  x  2  x  0.66  |  °*U* 

y  2.80  (-0.40  +  2.40  x  3)J 


For  comparison,  air  in  both  ohambors  give  w—  =  11.7  for  y  —  3 . 

e o 
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VI.  adiabatic  theory  of  the  shock  tube 

An  explicit  expression  for  P  oan  be  obtained  fran  an  adiabatio 
solution  of  the  bursting  diaphragm  problem  whioh  is  approximately 
correct  and  is  convenient  for  calculational  purposes.  The  errors  in 
the  approximation  are  indicated  in  tabular  form  at  the  end  of  thiB 
section. 


The  line  of  reasoning  is  the  same  in  this  case  as  in  the  previous 
more  exact  solution  but  the  expression  for  the  par.tiole  velocity  be¬ 
hind  an  adiabatic  compression' Wave  is  substituted  for  that  behind  a 
shook  wave. 

We  have  from  equation  (l.8b)  that  the  particle  velocity  behind  a 
oompressional  wave  of  pressure  P  is t 


r  ZLiI 

2a  -  -  2  r 

u  =7^1  V")  “  1 


partiole  velocity  be-  (VI.l) 
hind  compression  wave 


and  from  equation  (1.9)  the  particle  velooity  behind  a  rarefaotion  wave 
advancing  into  a  region  of  pressure  P  and  velooity  of  sound  a.  isi 

■*.  r  .  ^i° 


1  - 


(VI. 2) 


Assume  a.  =  a  as  before  and  equate  the  particle  velocities  giving 

7^-1  y-r 

n  T r 


-  <|-> 


(VI.3) 


*  v  v 

Y~  °  p"  =yF” 


,1  -  (y^) 


=  y 


(VI. 3a) 


p  r-i  2  r 

-  2r,r-i 

(2  -  y  )  ' 

and  likewise  the  explicit  solution  for  y  is 

2  r 

...  M  g*  >A 


[p  ^  1 


(VI.4) 


(VI. 5) 


If  Y  ss  1.4  (air)  and  ire  let  ?C/?Q  s  w  (  then 


2  7  co 


^ ;  D7 


an  explioit  expression  for  y  \ 

in  terms  of  oo. 


„  _  .128  OJ 

y  pT^TTF 

V  < 

also  in  terms  of  the  same  variable 


(VI.  5b) 


The  results  of  equation  (V.4)  give 


(VI  .4a) 


r,  ty-iyrr 
L  /  7(1  +  6  y)J 


and  CO 


l/7  = 


7U+ Tyj 


Substituting  this  value  of  CO  into  equation  (VI.5a)  and  extracting  the 
root,  we  have 

_  1/7  _  2  y1/7 

y‘  7 

/7(1  +  6y) 


so  that 


Ll  + 


3EI 

*6y)J 


Ratio  of  amplitudes  oal- 
-  oulated  for  the  compress-  /._  g\ 
'  ional  wave  under  conditions '  *  ' 

of  adiabatic  and  shook  wave 
formation. 


These  are  tabulated  as  functions  of  y 


(ya  -  y) 


X  100 


1  1.00  0 

2  1.01  1 

3  1.02  2 

4  1.033  3.3 

5  1.05  5.0 

6  1.07  7.0 

7  1.09  9.0 

8  1.115  11.5 

9  v  1.14  14.0 

10  1.17  17.0 

This  table  gives  an  idea  of  the  range  of  shook  pressures  over  which 
the  adiabatic  expression  for  y  may  be  used  without  too  great  an  error  in 
results. 
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VII.  THE  PRODUCTION  OF  ARTIFICIAL  BLAST  WAVES 

(The  calculation  of  the  relative  dimensions  of  the  tube  to  pro¬ 
duce  shock  waves  with  a  peaked  top.) 

The  shape  of  the  shook  Waves  produced  by  the  bursting  of  a  dia¬ 
phragm  in  a  tube  will  ohange  as  the  wave  progresses  down  the  tube  as 
shown  in  the  following  series  of  sketches. 


t«  *  ■  !bf  .  i 

t«  /  -«*y>  - 


Examples i  y  =  2 

ff  ~  11.3° 

fl'  «  -  14.3° 

angle  of  P  =  PQ 
in  the  trough. 


It  is  obvious  that  at  some  point  down  the  tube  the  rarefaotion 
tip  will  catch  up  with  the  shock  front  and  a  peaked  wave  will  be  pro¬ 
duced.  It  is  the  position  of  this  point  which  will  be  sought  in  the 
analysis.  The  sequence  of  events  in  the  tube  after  the  diaphragm  has 
broken  consists  of  the  production  of  a  shook  wave  traveling  down  the 
expansion  tube  with  a  velocity  U,  followed  by  a  temperature  dis¬ 
continuity  traveling  at  a  velooity  u,  and  a  rarefaction  wave  traveling 
back  into  the  compression  chamber.  The  initial  part  or  tip  of  the 
rarefaotion  wave  travels  back  with  the  velocity  of  sound  b  in  the 
compression  ohamber  which  is  assumed  to  be  equal  to.  aQ.  T?iis  tip  is  then 
reflected  from  the  closed  end  of  the  compression  chamber  and  travels 
for  a  short  distance  through  a  variable  density  region  until  it  reaohes 


the  gas  of  constant  pressure  P,  It  then  travels  down  the  compression 
chamber  and  expansion  tube  with  a  oonstant  velocity  b«  until  it 
reaches  the  boundary  between  the  cool  gas  that  was  originally  in  the 
compression  chamber  and  the  gas  passed  over  by  the  shock  wave  whioh 
has  been  compressed  and  heated.  This  boundary  travels  with  the 
particle  velooity  u  behind  the  shook  front.  When  the  rarefaction  tip 
reaches  this  boundary  it  proceeds  with  a  slightly  greater  velooity  bg 
until  it  eventually  catohes  up  with  the  shook  front. 

The  various  velocities  with  wtiich  the  different  sections  of  the 
shook  wave  travel  in  the  tube  are  derived  below i 

b0  =  S  (71 1.1) 

where  bQ  is  the  velocity  of  the  tip  into  the  compression  chamber. 
After  reflection  from  the  olosed  end  of  the  compression  ohaober  and 
after  reaohing  the  constant  pressure  region  P  the  velooity  of  the 
rarefaction  tip  is 


where 


and 


b^  =  u  +  (VII.  2) 

5a  (y  -  1) 

u  =  -  ■  °  . .  Equation  (III. 9a) 

J  7(1  +  6  y) 

1/7 

°1  "  *o  =  aQ  (i  “  .  (y  "  AL)  Equation  (V.4) 

1  °  Po  °  /7(1  +  6  y) 


so  that 

^  _  ,,  .  4(y  -  l)  velooity  of  tip  in  (VII.2a) 


The  velooity  of  the  tip  after  passing  the  gas  boundary  is 
bg  =  u  +  o 


(VII. 3) 


where  u  is  the  same  as  above  and 

Equation  (lll.8a)~-  .< _ __ 

so  that  _ 

V  _  .  f~5(y  -  D  +/7y(6  *  y)l  velocity  of  tip  (TTT.S^ 

2  °  L  v^ni  +  6  y)  J  “  hot 

After  breaking  the  diaphragm  the  velooity  of  the  trough  of  the 
rarefaotion  wave  baok  into  the  compression  chamber  will  be  vq  whioh  is 
the  looal  velooity  of  sound  minus  the  partiole  velooity.  So 

v  =  e.  -  u  (VII.4) 

O  X 
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so  that 


[*  w  _  . \  — |  volooity  of  trough  toward 

1  -  — -  W  -■  >-  the  olosed  end  of  the  ooia-  (VXI.4a) 
v/  7(l  +  6  y)  I  press  ion  chamber. 

It  may  be  noted  that  vn  =  0  when  y  =  2, 38  which  means  that  the 
trough  is-  stationary  at  this  shook  pressure.  It  moves  toward  the 
olosed  end  if  the  pressure  is  less  than  2.88  Pq  and  moves  toward  the  shook 
front  if  the  shook  pressure  exceeds  this  value; 


After  reflection  and  after  the  tip  of  the  rarefaction  wave  has 
moved  out  of  the  variable  density  region,  the  trough  of  the  rare- 
faotion  wave  has  a  velooity  v.  whioh  is  the  local  velocity  of  Bound  in 
a  gas  oooled  by  adiabatio'  expansion  from  a  pressure  P  to  a  pressure 

So 

r 


where 


but 


and 


so 


but 


Y1  =  °2 


(VII.5) 


pr 

ao  V5 
o 


1/7 


Pr 

(jr)  * 

■V  * 

■  2(y-)  - 

V 

and 

0 

,2  v1/7 

°2  *  ao 

,  y  r,  „  1) 

l/7 

6J  '  «  - 

yXA 

- *7  rr~ — rr 

p  XA  p 

t r  \  .  t  c 


1/7 


P  _ 

ro  ~ 


1  - 


/ 7(1  4  6  y) 


so  that 


r.  =  a  fl  . 

1  0  L  ^7(1  +  6  y)  J 


from  equation  (V.4) 


velooity  of  the 

trough  in  cool  (VII. 6a) 

gas. 


After  crossing  the  gas  boundary  the  velooity  of  the  trough  is 
equal  to  the  velocity  of  sound  in  a  gas  cooled  from  a  pressure  P  and 
sound  velocity  c  to  a  pressure  P^,,  and  a  sound  velooity  a_  plus  the  .«■ 
residual  particle  velooity  which  is  no  longer  zero  after  orossing  the 
gas  boundary. 

Rrom  section  II  we  have  the  equations  which  give  the  value  of  Prg 
after  passing  through  the  gas  boundary  as 

P  ,  *A  2a  P.  1/7  0  -  a. 

<-£>  '  o hq  &  *  r-rq 
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4 


*1  c  ao 


P  i  A 

(p~) 

0 


,p  .V7  1/T 

and  from  adiabatic  theory  of  the  tube  we  have  fa—)  *=  2  -  y  ' 

re  ' 

eo  *1  B  ‘c^  "  y^> 

and  o  *=  aQ  ^  from  e9ua’fc^an  (HI.8a) 

IT  an  the  theory  of  reflection  of  adiabatio  waves  (seotion  XI)  we  have 

Prl  P0  .  2  1/7  P0  1//?  . 

hr)  =  2  -  (y^)  *=  ^ -  ^7-  since  (y-)  « 


from  equation  (III. 8a) 


Substituting  these  values  we  have  then 

Pr2  ^o  1/7  °  “ 

<-r>  “tttt  <3-2y1/7)+ — 


o  —  a. 


(vxx.r) 


So  the  velocity  of  the  trough  will  then  be 


where 


v2  =  u  “  U1  4  a3 


6a  (y  -  1) 


7(1  +  6  y)  y  7  y(6  +  y) 


u.  s=.5o  X  - 
A  1— 


pr2  lA 
(-£> 


the  partiole  velocity  behind 
the  shook  wave 


the  partiole  velocity  caused 
by  the  rarefaction  wave 


and  *5  =  c  (-|^) 

Substituting  values  for  these  quantities  we  then  have  that  the 
velooity  Vg  of  the  trough  in  the  hot  gas  is 

r  ..  _  ,  *  2a  (3  -  2  y1/7)  +  o  -  a,  V 


=  So{s?W? 


-  1  4 


2a  (3  •  2  yV')  +  o  -  a.  ^ 

<-  —  o— - Mjmi-e) 


The  position  of  the  shock  wave  in  the  tube  at  any  time  may  be 
graphically  portrayed  in  a  ohart  which  shows  the  locus  of  the  various 
points  as  a  funotion  of  time  and  distance  down  the  tube.  If  distance 
along  the  tube  from  the  diaphragm  divided  by  the  length  of  the  com¬ 
pression  ohamber  is  plotted  as  x  along  the  absoissa  and  if  the 
quantity  a  t/LQ  is  plotted  along  the  ordinate,  then  straight  lines  may 

be  drawn  (representing  the  locus  of  the  points  as  a  funotion  of  time) 
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with  the  position  of  the  diaphragm  as  the  origin  and  at  angles  0.,  0g, 
ete.  with  the  absoissa.  The  tangents  of  the  angles  0^  0g,  etc,,  are 
the  inverse  functions  of  the  velocities  of  the  pointer  We  then  have 
defined  for  the  various  components  representing  shook  Velocity,  tr  ough 
velooity,  rarefaction  tip  velooity,  etc,,  the  tangents  of  the  angles 


of  their  looi 

as 

follows* 

a 

tan 

*1 

0 

u 

shock  front 

(VII.9) 

tan 

®2 

- 

a  ;*■ 

o  '  • 

u 

gas  boundary 

(VTI. 10) 

tan 

®3 

«K 

*D 

u 

rarefaction  tip  in 

(VII. 11) 

0 

0 

0 

ii 

compression  chamber 

tan 

«4 

e 

tlo 

TJT 

rarefaotion  tip  in  oool 

(VII. 12) 

tan 

*5 

s= 

b2 

rarefaction  tip  in  hot  gas 

(VI I. 13) 

tan 

e6 

“5 

a 

_£L 

trough  into  compression 

(VII. 14) 

¥ 

0 

chamber 

tan 

«7 

zz 

±2 

V1 

trough  in  oool  gas 

(VI 1. 15) 

tan 

08 

= 

1q 

v„ 

trough  in  hot  gas 

(VI  1.16) 

A  time  distance  chart  for  a  shock  wave  having  a  pressure  ratio 
y  =  2  is  shown  on  a  following  page,  (Graph  7.)  From  it  a  groat  deal 
of  information  may  be  gleaned.  For  example,  the  length  of  the  ex¬ 
pansion  tube  necessary  to  insure  that  the  shock  wave  have  a  pointed 
top  like  a  blast  wave  can  be  found  at  the  intersection  of  the  locus  of 
the  shock  front  (line  U)  and  the  locus  of  the  rarefaction  tip  (line  b_). 

In  this  case  it  appears  that  the  length  of  the  expansion  chamber  should  be 
15  times  the  length  of  the  compression  chamber  to  obtain  this  condition. 
The  duration  of  the  shock  wave  from  front  to  trough  can  be  found  in 
terms  of  aQt/Lc  at  any  distance  x  along  the  tube  by  measuring  the 

vertical  distance  from  line  U  to  the  line  v„.  In  similar  fashion  the 
duration  of  the  flat  top  of  the  shock  wave  at  any  point  x  is  found  by 
a  measurement  of  the  vertioal  distanoe  from  line  U  to  line  b„. 

Another  method  of  showing  the  ohange  in  shape  of  the  shock  wave 
while  passing  through  a  homogeneous  medium  is  to  show  the  looi  of  the 
points  relative  to  the  shock  front  as  in  the  first  page  of  this  seotion. 

In  this  case  the  abscissa  is  the  length-  or  duration  of  the  shook  wave 
while  the  corresponding  ordinate  is  the  distance  along  the  tube  or 
time.  In  this  caso  the  angle  of  tho  locus  with  the  ordinate  is  <X  where 

tan  d  =  J  -  1  (VII. 17) 

where  v  is  the  velooity  of  the  component  under  consideration  and  U  is 
the  shook  wave  velocity. 
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II 


HiiHgigfl 

SB 


The  oritioal  length  of  the  expansion  tube  whioh  first  allows  the 
rarefaction  wave  to  reaoh  the  shook  wave  may  be  solved  by  calculation 
rather  than  graphically  if  we  proceed  as  follows i 


After  the  diaphragm  is  broken  the  boundary  between  oool  and  hot 
moves  down  the  expansion  chamber  with  a  ■yelooity  u.  where 

5a  (y  -  1) 

-  -===== 

✓  7(1  6  y) 

The  distanoe  that  it  has  moved  at  the  instant  that  rarefaotion 
break  reaohea  the  end  of  the  compression  chamber  is 

L0 

d  su  l  but  T  »  — 


(vii.ie) 


so 


(VII. 19) 


The  length  of  the  oool  gas  oolumn  at  this  instant  is 


n  = 


Lo  + 


Ln  (1  4  ■  ■  pjy  -JL)  (VII, 20) 
/7(1  +  6  yj 


The  time  of  travel  neoe^sary  for  the  rarefaction  break  to  reaoh  the  oool 
gas  boundary  is 


s(y  -  l) 


1  +  6  y) 


y  -  i) 


/ 7(1  4  6  y)  j 


(VII. 21) 


The  distanoe  x  from  the  diaphragm  at  whioh  the  boundary  is  reached  is 

uL  L 

x  =  u  T  4  — =  u  (T  4  —.)  (VII. 22) 

cL  R  _ 

O  O 

Substituting  from  equation  (VII.21)  we  have 

x  =  6,L°(y : 1}  Hi*  i)  (vn.23) 

✓  7(1  4  6  y)  ^7(1  4  6  y)  -  (y  -  l) 

The  length  of  the  hot  gas  column  at  the  time  the  cool  gas  boundary  has 
progressed  a  distanoe  x  down  the  tube  is 

%s=  t(U  -  u)  where  T  =  ^ 


* 
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5  =  x(-  -  l) 


6  +  y  H 


(VI I. 24) 


The  additional  distance  that  the  r&refaotion  break  must  traverse  to 
pass  through  the  hot  gas  and  reaoh  the  shook  front  is 


*i  ■ 


where 


*1  *  ttt* 

1  S 


1  +  6 


2  6(y  -  1)  +  /  7  y(6  +  y) 


SO  ®  X 


fe  +  y  Cs(y  -  1)  +  >/7  y(6  +y>;  J  1 
l  5(y  -  1)  £V 7  y(6  +  y)  -  (6  +  y)]  J 


(VXI.  26) 


The  total  distanoe  from  the  diaphragm  at  whioh  the  rarefaotian  break 
overtakes  the  shook  wave  is  the  sum  of  these  two  distanoes  x  and  x^ 

.0  . ,  -  -  |~ '  y)  *  6(y  -  1)' 

1  (_  My  1)  V7y(6  +  7).  (6  +  y). 

We  then  set  ff  = 


£  *=  1  + 


r7  y(6  +  y)  +  5(y  -  l] 


so  that 


V  " 


5(y  -  l)[y  7  y(6  +  y^  -  (6  +  y)] 


for  the  oritioal  length  of  tube* 


The  values  are  tabulated  below. 


y 

V*o 

P0S  lb 

? 

< 

L* 

1.25 

1.558 

8.2 

0.3566 

75,85 

27.0 

1.50 

2.29 

18.97 

0.712 

25.81 

16.4 

1.75 

3.20 

32.35 

1.064 

15.0 

15.95 

2.00 

4.36 

49.4 

1.418 

10.64 

15.1 

2,50 

7.28 

92,3 

2.12 

7.03 

14.9 

3.00 

11.40 

153.0 

2.82 

5.50 

15.5 

4.00 

24.16 

340.0 

4.255 

4.142 

17.60 

These  oritioal  lengths  in  termB  of  y  and  PQa  are  plotted  on  -the  next 
pages*  (Graphs  8  and  9.) 
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Graph  8 


Critioal  Batio  of  Lengths  of  Expansion  Chamber  to  Length  of  Compression 
Chamber  as  a  Function  of  Shook  Strength  (y) 
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VIII 

THE  IMPULSE  IN  THE  SHOCK  WAVE  AT  THE  POINT  OP  CATCH-UP  OF  THE 

RAREFACTION  WAVE 


A  solution  of  this  problem  may  be  attempted  by  assuming  that  the 
energy  released  from  the  compression  chamber  by  the  expansion  of  the 
gases  from  a  pressure  P  to  a  pressure  P  ,  the  pressure  in  the  trough 
of  the  reflected  rarefaStion,  is  essentially  conserved  and  is  the 
energy  available  to  do  work  in  the  gas  behind  the  shock  wave.  The 
energy  in  the  shock  wave  is  conceived  to  be  the  internal  energy  of  the 
compressed  gases  in  the  shook  wave  if  these  gases  are  expanded  adia- 
batically  to  the  pressure  Pr  of  the  rarefaction  wave  following  the 
shook  wave  plus  the  kinetic  energy  of  the  gas  partioles  in  the  shook 
wave.  The  assumption  of  conservation  of  energy  will  be  examined  later 
to  determine  the  error  introduced  by  it. 

If  the  conoept  of  the  energy  in  the  shock  wave  outlined  above  is 
adopted  it  will  then  be  permissible  to  assume  &  relationship  between 
the  energy  in  the  shock  wave  and  the  impulse  associated  with  it  of  the 
following  form 

(VIII. 1) 

where  E  «b  energy  of  the  shock  wava 

I  =  impulse  of  the  shock  wave  as  measured  by  a  gauge  side  on 
U  *=  velocity  of  the  shock  wave 

and  k  =  a  dimensionless  factor  to  be  determined  later.  It  is 
a  function  of  the  shock  strength  and  the  shape  of  the  pressure-time 
ourve  of  tha  Bhook  wave. 

Equation  (VIII.l)  will  be  used  to  calculate  the  impulse  of  the 
shock  wave  when  we  have  evaluated  k  and  E,  the  energy  in  the  shock. 


The  work  done  by  a  unit  mass  of  gas  in  adiabatioally  expanding  from 
a  pressure  Pq  to  a  pressure  Pp  is 


but 


so 


e  -  e 

o  r 


/°o  _ 

7T  -<r,) 

~  r  r 


g  -  e 

o  r 


-  t*  [k  ■  k]  -  r*  k  [*  - % 

i/y 


1  -c 

r^r  *r 


r 

p  r 

r  1  -<r>  J 

c  L  .  r  — 1 


energy  per  unit  (VIII.2a) 
mass  of  gas. 


The  available  energy  per  unit  volume  of  gas  at  a  density  /°  in  a  tube 
of  length  Lq  ana  unit  area  is  G 

e  =  /°  L  (e  -  e  ) 

'  o  O  '  C  Tr 
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so 


0  = 


1 


and  if  r  «=  1.4 


PL 

O  0 

y=x 


I  PoLc 


(VIII.3) 


This  is  the  energy  available  for  adiabatio  expansion  in  a  cylinder  of 
gas  of  unit  cross  section,  of  length  LQ,  at  a  pressure  PQ  and  expanded 
to  a  pressure  ?r. 

If  P  is  the  peak  pressure  in  the  shook  'wave  and  P  is  the  initial 
pressure  in  the  compression  and  expansion  ohamber  we  tnen  define  the 
ratio 

P/P  =  y  and  P_/P  =  O) 

*0  o#  o 


If  we  consider  the  initial  formation  of  the  wave  to  be  adiabatio 
we  oan  get  an  explioit  expression  for  y  in  terms  of  60  whioh  is i 


zgL l 

W1/7  +  i 


equation  (VI.5a)  (VIII.4) 


This  assumption  whioh  is  made  for  computational  convenience  results 
in  an  error  which  may  be  determined  for  any  range  of  shook  pressures 
from  the  table  at  the  end  of  section  VI.  For  the  range  of  pressures 
oonsidered  here  it  may  be  responsible,  for  a  maximum  error  of  about  2 
percent  at  the  highest  shook  pressure. 


Team  the  theory  of  reflection  of  adiabatio  waves  we  find  the  ratio 

p y 


=  2  y1/-  -  to1/7  equation  (11.8b)  (VIII.5) 


consequently  in  terms  of  CO 

p 

$ 


but 


i 

^  *  i 


* 


i A  P  iA  P  iA 

=  (/)  tf) 

o  r 


(VI II.  6) 


so 


r 


-*A  1/7  ; 

-  (s \  %  x  ) 

+  1 


s  —  p  L 
2  c 


8(<0 


1/7  . 


(VI II.  7) 


then 


E 


(rai.  e) 


(VIII. 8a) 


which  is  also  equal  to 


[■ 


<o(tO 


lA . 


1) 


E=20*°l°  IT^TT? 


] 


This  equation  (VTII.8a)  then  gives  the  value  in  terms  of  6>of  the 
available  energy  in  a  tube  of  unit  cross  section  and  length  Lq  when 
expanded  to  a  pressure  P  whioh  is  determined  by  the  reflection  of  the 
rarefaction  wave  produoeS  by  a  bursting  diaphragm. 

S  < 

If  a  shook  wave  of  pressure  P  is  traveling  into  undisturbed  air 
of  pressure  P  the  velooity  of  the  shook  wave  is  tJ,  where  U  is  given  by 
the  equation 


../i 


4  6  y 


(VIII, 9) 


Again  using  the  explicit  expression  for  y  in  terms  of  co  given 
above  for  the  adiabatic  case  we  have 


U  =  ■ 


r?  L1  +<  j/7 + D7 


768  60 


■] 


1/2 


velocity  of  shook 

wave  in  terms  of 

tO  K  P  /P 
o'  o 


(VI II. 10) 


Vfe  now  have  expressions  for  two  of  the  faotors  involved  in  equation 
(VIII.l)  and  it  remains  to  evaluate  the  factor  k  in  order  to  arrive  at 
an  answer.  One  might  infer  from  an  analogy  with  the  mechanical  case 
where  k  =  2  that  the  factor  might  lie  in  the  neighborhood  of  this  value. 
The  oonoept  of  stored  energy  in  the  shook  wave  coupled  with  some  reason** 
able  assumptions  as  to  the  shape  of  the  pressure-time  ourve  of  the  shook 
wave  enables  U3  to  make  a  reasonably  accurate  determination  of  this 
faotor  which  does  turn  out  to  be  of  the  order  of  magnitude  of  2. 

We  proceed  to  evaluate  k  in  the  following  manner. 

From  equation  (VIII.l)  we  see  that 

k  =  ™  for  any  blast  wave.  (VIII. ll) 


The  energy  in  a  shock  wave  consists  partly  of  potential  and  partly 
of  kinetic  energy,,  the  proportions  of  which  change  with  the  amplitude 
of  the  shock.  The  potential  energy  is  considered  to  be  the  available 
pressure  energy  of  the  gas  in  the  shock .wave,  while  the  kinetic  energy 
is  the  energy  of  motion  of  the  particles  in  the  shock  wave. 

The  available  pressure  energy  of  a  slice  of  gas  Bomewhere  behind 
the  shock  front  is 

ABi  s  -  7^)  =  “7T-  1)  ( viii. 12) 
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but  P/fcQ  e  y  and  /°//°0  -  x 


1 

y  -1 


-  1) 


Energy  per  unit  mass  of  gas 
if  the  gas  were  expanded 
adiabatically  to  a  pressure 
PQ  and  density  /°  * 


(7111.12a) 


The  energy  per  unit  volume  is  /&  &  Ej  and  is 


^Bpc  p=T  *0(y-*)  (7111*18) 

and  if  X  **  1*4  (air) 

then  ABp  c  |  P0  ^  (Till. 13a) 

From  equation  (III. 5)  seoticn  III,  we  find  that  the  ratio  of 
densities  x  in  a  shook  wave  is 


x 


14  6 
6  + 


J l 
y 


5 

2 


p  (2 _ Z_l) 

ov  6  4  y 


Potential  energy  in  a  slice  of 
gas  of  unit  oross  section  and  /___ 
unit  length  at  pressure  P  and  '  •  5  / 

density  /°  • 


The  total  potential  energy  in  the  blast  wave  is  thon 


E 

P 


dx 


(VII I. 14) 


The  kinetio  energy  per  slioe  of  length  dx  and  unit  area  is 


AE2  =  4  P* 


where  u  is  the  partiole  velocity  .R\ 

and  /o  is  the  density  VVAAJ,*A0' 


(7111.16a) 


but 


so 


2 


u 


26  a?'(y  -  l)2 
7(1+  6  y)“ 


1  -  26^  (y  -  l).2,  d* 

2  f  o  7  6  +  y 


1  +  6  y 
6  +  y 


( VIII. 16b) 
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And  the  total  kinetic  energy  in  the  wave  is  then 


,  = 28 

*k  I?  V  6  +  j 


dx 


(VIII. 16) 


hut  a2  c  1  -g- 

5  So 


BO 


i  p 


(vinaea) 


6  •?  f*  „(y,  -  U2_i* 

O J  6  +  y 

•/0 

The  total  available  energy  in  the  shock  wave  is  then  the  sum  of  the 
kinetio  and  potential  energies  which  is 


b.  «  /  £. i A..4  .(y  --1)2 **. 

t  2  ro  /  6  4  y 

'o 


*=  5  P 


// 

rfy-lL?*  total  energy  in  the  (VIII,17) 

6  4  y  shook  wave 


Now  if  the  shook  wave  changes  in  duration  slowly  enough  while 
passing  over  its  own  length  we  may  put 


dx  =  U  dt 


(VIII. 18) 


ao  that  these  equations  may  be  written  as 


6  4  y 


(y  -  D2 

6  4  y 

xk  -_U, 

6  4  y  . 


dt 


dt 


dt 


(VIII, 19) 

(VIII. 20) 

(VIII,  21) 


We  may  substitute  z  =  y  -  1  as  a  variable  and  make  an  assumption 
as  to  the  rate  of  ohange  of  z  with  time.  Oscillograms  of  the  pressure- 
time  variation  show  that  the  wave  shape  approaches  triangular  form  at 
the  lower  pressures  while  at  the  higher  pressures  it  approaches  an 
exponential  form. 
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Assume  that  the  nave  shape  is  exponential  and  that 

-  -  1  rtr. 

*  =  *  e  then  dt  =  -  £  2L 

when  e  =  0,  t  -  o©,  z  =  z ,  t  *=  0 

so  that  equation  (VIII.21)  may  be  written  after  substituting  for  y  and 
dt  Jz 

^  PoU  f  (s  +  l)  dz  , _ _ * 

Bt - WT“  J  A  7  V  z -  (VIII. 22) 

.  .  5Po"  C  .4.  „  f  « 

Kt  -  -TT-  J  m  m 

5  p  d  *:  s 

*=  |7  +  *  -  7  log  (7  +  z/j  +jlog  (7  4  t)  j 

5  P  U  T  -  “I 

=  [S  .  6  log  (1  4  £)J  (VIII. 23) 


5  PD 
o 


The  impulse  measured  side  on  in  a  blast  wave  is  defined  as 


-  /».«  -  *0  (  (y  -  1) 

Jrs  Vi 


dt  c  p 


but  dt  «  V 


I  = 


dz  - 


(VIII. 23) 


*  dt  (VIII. 24) 


(VIII. 25) 


IU  = 


o  ,  IU 

and  k  =  -jr 


therefore  the  value  of  the  faotor  k  may  be  written  as 


k  = 


5  \  z  -  6  log  e  (14  iM 


(VIII. 26) 


The  term  log  (1  4  y)  may  be  expanded  in  a  series  to  give  the 
approximate  formula 


c /,  ,  3z  2z  3z  % 

5(1  4  T  “  T?  +  656  ••••> 


(VIII. 26a) 


which  has  the  value  7/5  when  z  — -^0.  This  is  of  the  order  of  magnitude 
of  the  meohanioal  factor  2  as  mentioned  before. 
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k  may 


X  2  B  c o 

If  the  explioit  relationship  y  ~  y 

■  be  written  (co1' '+  1/ 


is  used  the  faotor 


128  Co  -  (to1/7*  l)7 


6  £128  «  l)7{l  +  8  los8  (f  *  j 


(VIII,26b) 


Recapitulating  the  formulas  wd'have 


IcT 


128  -  (to1/7  4  l)7 


S  J*128  <*>  -(^A*.  l)V{l  4  6  log^  [j  4 


a  r 

u  =  -S  14 

/7  L 


- _ _  1//?  -  1) 

B  -  20  poLo  I  /rx/r—^z 


(a>v '  +  i)‘ 


7(tx) 


12§'co  Tl 

J 


(VIII.  n) 
(VIII. 26b) 


768  60 


(fiA'  4  1) 


7 


1/2 


(VIII.8a) 


(viii  ao) 


Combining  these  factors  we  then  have  the  expression  for  the  impulse 
in  the  shock  wave  in  terms  of  the  dimensionless  quantity  Ia^/P^l^. 

This  equation  will  hold  at  the  point  of  catoh  up  of  the  rarefaotion  wave 
with  the  shook  front  if  the  energy  lost  by  the  shock  wave  in  passing 
down  the  oritical  length  of  tube  is  ignored.  This  will  result  in  an 
error  of  less  than  4  percent  in  the  calculated  value  of  the  impulse. 

We  have  then 


Iao 

PoLc 


[l28  60  -  (CO1/7  4  l)7, 


(VIII.27) 


||M— <«VW{»«  1.*. 


r  768  co  -|1  /i 

L1  +  T^7J  ' 


This  relation  is  plotted  in  terms  of  (  o>  -  l)#  Graph  10. 


IX.  THE  CHANGE  OF  ENERGY  AS  THE  SHOCK  WAVE  PROGRESSES  ALONG 

THE  EXPANSION  CHAMBER 

When  a  shock  wave  passes  over  a  mass  of  gas  it  changes  the  entropy 
of  the  gas  and  leaves  it  in  a  different  energy  state  in  general  than  it  was 
originally.  This  ohange  in  energy  can  be  found  by  differentiating  with 
respeot  to  time  the  equations  which  represent  the  energy  of  the  shook 

wave. 
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4VW 


Graph  10 


Impulse  in  Shook  Wave  as  a  Function  of  Exoess  Pressure  in  Compression 

Chamber 


I  *=  Impulse  in  lb/seo/aq  in. 

«=  Velooity  of  sound  ft/seo* 

P  -  Pressure  in  expansion  ohamber  pounds  (atmospherio  normally) 
Lq  «=  Length  of  expansion  ohamber  feet 

w  *=  p  /p  y 

o'  0  / 

P  =  Pressure  in  compression  ohamber  lb/sq  in.  (total)  S' 
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If  one  considers  a  shook  wave  in  which  the  pressures,  densities,  sad 
partiole  velocities  have  values  indicated  by  the  subscripts  on  the  dia¬ 
gram  as  follows. 


P_  /Ort  u  =  0 
0^0  o 


the  coordinates  at  any  instant  of  time  of  the  boundaries  between  these 
regions  have  the  values  XQ,  and  Xg.  The  total  energy  in  the  shook 

wave  then  is  given  by  the  integral  of  the  bubi  of  the  kinetio  energies  of 
the  partioles  and  the  inorease  of  internal  energy  of  the  gas.  This  is 
given  below  as 

-X, 

I  A 

(IX.1) 

a  v 

'2 

The  change  of  energy  with  respect  to  time  is  the  derivative  of  this 
expression  with  respect  to  time  and  is 


X  /-X, 

B  *  f  A  (i  U,2+  Er  E0)dx  4/  uV  E  -  B0)dx 
•4,  JxP 


3E  “i2+ 


*/ 


1  “o'3t 


3xo  z1  2 ,  r  -  X^x2 

"  ^2%  u2  E2‘Eo%? 


3 

at 


u2  +  E  -  E  )dx 


(n.2) 


Now  — -r  —  U  the  velocity  of  the  shook  wave 

a  t 
and 


ax„ 


*  +  U„ 


~  ~2  '  “2 


where  a,,  is  the  velocity  of  sound  in 

the  trough  and  u,  is  the  partiole  (IX«3) 

velocity  in  the  trough. 


So  the  rate  of  energy  change  is 

■ff  u12+E1‘B<.>-  u22+  V  V 


+  /  /°(|  U2  +  E  -  Bo)  dx 


(IX.4) 
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The  last  term  in  the  integral  sign  is  equal  to 

.f  -ft (?  +lp  ’  I  jir/a'>  **  * 

xz 

f </”*  It.  +iuZ  -14 + 1  H  * !  -¥?)  «*  <“.«•> 

✓y  O 

X2 

t  < 

but  from  the  equation  of  motion  and  continuity  we  hare 


■m 


-yOSS.  U 

rax  ax 


(IX, 5) 


9i  ^  $u  1  3P 

at  "u  a*  “  ax’ 


(IX.6) 


Substituting  uheso  values  the  integral  becomes 


(„  3  xju2  —  -  iu3  -IS£  +  £  I°_  ...  5  Po  -  .&P 

'  3T^U  5x  2  ^§>x  2  pf  ax  2/®0  Sx 


/ 


-U*JLu-2£  -  5  _ap  _  9u,  .  r*  a>  ,1  - 3\ 

2  a/)u  ax  2  3/0  ^  ax'  ^  y  9x  'T  ^ 


p 

+  -2.(6 
ax  '2  /« 


/>u)  -  7  (  9P  5  _SP  .CM*) 
"  2  '  ax  7  ®P  r  ■a’r/ 


3U| 

9xJ 


(IX*  7) 


.  ,  5  ap  ~  _  5  7  P  ^,-p 

but  7“a^/^7*5’;o*/0=-p 


(IX. 8) 


So  thon  performing  the  integration  and  substituting  in  the  limits  we  find 
that  it  is  equal  to 

1  jQ  3  5  o  jn  7  D  ,  1  jQ  3  5  Po  /j 

**  ■g’  r  \  +  7T  rg—  r\  U1  **  ‘1  U1  "*  2  '  Z  u2  "If  !3“  '  2  w2 

’  o  *  o 

+  |  ?2  u2  (IX.9) 


Then  adding  all  terms  we  have 


9E  _ 

at 


=  /°1»4  “l2+  V  Eo>  -  firth  V  Eo>  -  «r+  V  V 


rl  ..  2 


1  2, 


-  PlWl  +  P2u2 


(IX, 10) 
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and  since  /°  .(TJ  -  u.)  *  A  V 
then  * 


by  the  equation  for  conservation  of  mass. 


ff  =  /°  oD(I  V  Eo>  -  Pl“l  -  /V2<1  “22+  V  V  +  V2  <«•«> 

However  as  a  consequence  of  the  conservation  of  energy  as  defined 
in  the  Rankine-Hugoniot  equations  the  first  two  terms  oanoel  so  we  are 
left  with  the  expression  for  energy  change  as  follows* 


9  B 

a* 


-  A 


2*2 


,1  2 

%  u2 


V 


Eo> 


P2U2 


where  Ug  must  be  given 

the  proper  veotor  sign.  '  *  ' 


This  equation  expresses  the  idea  that  the  energy  change  as  a 
funotion  of  time  is  the  difference  between  the  total  energies  of  the 
gas  as  it  enters  and  leaves  the  shook  wave  plus  the  work  done  on  or  by 
the  gas  after  it  leaves  the  shook  wave  depending  on  whether  ug  has  a 
positive  or  negative  sign. 


Let  us  consider  the  shook  wave  after  the  rarefaotion  trough  has 
passed  the  gas  boundary.  The  reason  for  doing  so  is  that  until  this 
stage  is  reached  the  gas  in  the  compression  chamber  is  still  giving  up 
energy  to  the  shook  wave  but  after  this  stage  the  shock  wave  etarts.-to 
dissipate  and  to  feed  baok  a  small  amount  of  its  energy  into  the  eom- 
pression  chamber  in  the  form  of  a  compression  wavs. 


After  this  point 


u 


2  *  u 


where 


so 


u  =  5c 


u„  =  5o 


c  -  a. 


XT-1. 


✓  7  y(6  +  y ) 


y  -  i 


y 7  y(6  +  yj 


TEL ±  y). 

1  +  6  y 


and  \ir  ®  5c. 


1/7 


-  1  + 


1/7 


i  -  (4) 


where  y  =  — 
o 


(12.13) 

(IX.  M) 

(IX.15) 


Now  from  the  theory  of  reflection  and  transmission  of  rarefaction  waves 
as  given  in  seotion  II,  we  have 

■  p  1/7 

2a,  r.  c  »  a, 

hh  +— TT 


Pj> 

hi) 


1/7 


where  o  =  speed  of  sound  in  the  shock  wave 

a.  =  speed  of  sound  in  gas  that  has  been  oooled  by  expansion  from 
pressure  P„  to  P„.  • 

o  rl 


Tr am  seotian  II  vre  have 


r. 

(-pi) 


1/7  P  iA 

*  2  -  <f  > 


(a.  17) 


hut 


p  1/7 


iA 


i  -  _(y  -  I) 

</  7(1  +  6  y) 


(IX. 18) 


80 


(-tt) 


/7(1  4  6  y)  2(y  l) 
/'7(l  4  6  y)  -  (y  -  l) 
p  l/7  p 


(IX*  19) 


1/7 


and 


al  - 


(-r-)  -shr»  %> 

0 

[l.  *<*--  1)  1 

’  1_  /  7(1  +  6  y)  J 


(II.  20) 


The  ratio  of  densities  /°  /  y°g  oan  be  found  by  considering  the  gas 
to  be  compressed  aooording  to  tfie  Hankine-Eugoniot  equation  of  state  and 
to  expand  aooording  to  the  adiabatio  equation  of  state  so  that 


1  *  «  r 

6  4  y 


and 


/v 

~F 


p  pA 

<■#) 


then 


then 


a2  a. 

*  T 


•  c  t  o  _  /  2» 

-Tr  •  -zr  =  vy/ 
0 


5/7 


4 


5/7 


(IX.21) 


A0 


-5 A 


p  -5A- 

*  A 


then  /°2  (gr) 


5/7  *  V  (^)  =y5/7.(fr^) .  (p^)  (IX.22) 

P  5/7 

/V-^o  (>r>  •  VrS'  <LT*23> 


5/7  r 


O  ■ 


14  6  y 

Ly6^(6  +  y) 


The  ohange  of  internal  energy  Eg  -  Eq  is  equal  to 
_  «  5/P2  ?°l  _  5  Po  fP2  ^0 

*2  -  Bo  "*k  *  K  <7l  •  7®I " 17 


(IX.  24) 
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Substituting  for  f*  J  P  2  have 


h-*a 


6  Po 

7  ^ 


y6/7<6  +  v>  ,\Z/1 

nr+Ty  (-r) 


•*j 


1/7 


and 


ft2  5=0  <T>  =  % 


/y(6  +  y)  ,P2.  ^ 
y  l  +  o  y  (T> 


Example  y 


0 

a  /zLG.J.  y)  = 

°V  1  +  6  y 

ax 

«=  0.790  a 

0 

P  i  A 
<V> 

*  0.902 

p.  1/7 

(y) 

0 

«=  0.9965 

u 

—  -*0.0194  a^ 

'oVT3 '  "  1,11  V 


2 

1  2 
2  u2 


0.000263 


7r 

'  o 


/°2  *  0.978  /O' 


p2 

T 

0 


0.977 


(IX.25) 


(IX.26) 


E„  -  E 
2  o 


0.0025  ■£■ 


a*  —  1,00  a 


then 

9E 


3  b  r 

9t  =  “Vo  L0*978  *  1«°°  (0.000263+  0,0025)  -  (0.977  x  0.0194)]  B 

Vo  (°»0165) 

The  total  dissipation  Ed  =  t  =  0,0165  P^t 

Fran  the  diagram  of  section  VII  we  find  that  the  transit  time  of  the 
rarefaction  trough  from  the  gas  boundary  to  the  point  of  catch  up  is 
L 

11  ~  so  that  a  t  s  111  , 
a0  oo 


Therefore  the  total  dissipation  E,  ia  0.181  P  L  . 

d  oo 
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The  energy  originally  in  the  wave  is 


B  =  10  PL 


y(yl/V;  ll 

00  LiTTtj 


=  4.03  PL 
o  o 


(IX. 27) 


so  that  the  ratio  of  the  energy  dissipated  to  the  total  original  energy 
1b 


Kd  _  0.181 

t  c  T^r 


0.045  2  5% 


This  gives  a  measure  of  the  aoourooy  of  the  calculations  of  impulse 
■which  are  based  among  other  things  upon  the  assumption  of  conservation  of 
energy  of  the  shock  wave.  Presumably  then  at  the  pressure  level  y  e  2 
the  error  in  the  oaloulations  should  be  of  the  order  of  5  percent. 

X.  AN  EXPERIMENTAL  DETERMINATION  OP  THE  POINT  OP  CATCH  UP  OP 
THE  RAREFACTION  WAVE  AND  THE  IMPULSE  OF  THE  SHOCK  WAVE  IN  A  TUBE 


ABSTRACT 

Experiments  are  reported  here  which  show  that! 

(A)  The  impulse  in  the  shock  wave  in  a  tube  at  the  critioal 
distance  is  very  nearly  proportional  to  the  excess  chamber  pressure  and 
oan  be  caloulated  by  considerations  of  the  energy  in  the  shook  wave. 

(B)  The  length  of  expansion  tube  neoessary  to  allow  the  rare- 
faotion  wave  to  oatoh  the  shook  front  and  to  produoe  a  peaked  shock  wave 
oan  be  found  from  a  knowledge  of  the  various  wave  velocities  and  is  of 
the  order  of  magnitude  of  fifteen  times  the  compression  chamber  length 
at  the  higher  shook  pressures. 


Note!  The  complete  paper  on  this  subjeot  will  be  found  in  AES 
7  February  1945. 
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APPENDIX  A 


A  SIMPLE  DERIVATION  OF  THE  EQUATION  FOR  THE  TOTAL  ENERGY  IN  A  SHOCK  WAVE 


Assume  a  tube  of  unit  oross  section  open  to  the  atmospheric  pressure 

P  at  one  end  and  closed  at  tho  othor  end  by  a  weightless  piston.  The 

air  behind  the  piston  is  evacuated  to  avoid  any  complications  due  to  the 

formation  of  rarefaotion  waves  so  that  the  statio  pressure  on  tho  piston 

is  PQ,  Now  let  the  piston  be  aotod  upon  by  a  suddenly  applied  uniform 

force  P  which  causes  it  to  move  with  a  uniform  velocity  u  until  the 

piston  has  moved  a  distance  xq„  The  total  work  done  by  the  piston  is 

then  Px  . 

o 

The  motion  of  the  piston  will  produce  a  region  ox-  oompresBed  gas 
ahead  of  the  piston  whose  pressure  is  P  and  whose  boundary  will  be  ft 
shook  flrant  of  velocity  U.  The  velocity  of  the  piston  will  be  u,  the 
velooity  of  the  gas  particles  behind  a  shock  lYont  of  velooity  U  and 
pressure  P.  The  shook  front  will  then  move  a  distance  Ut  while  the 
piston  moves  a  distance  ut.  The  length  of  the  shook  wave  at  the  time  t 

Xo 

will  be  (U  -  u)t.  Let  t  «  ~  then  the  length  of  the  shook  wave 
at  the  time  the  piston  has  moved  a  distance  xq  will  be 


x 


1 


(A.l) 


The  energy  in  the  shock  wave  will  then  be  equal  to  the  energy  per 
unit  volume  times  the  volume  of  the  shook  wave. 

Since  the  tube  is  of  unit  cross  seotion  then  this  is 

<§>  •  *i  the  total  energy  in  the  shook  wave  (A. 2) 


This  in  turn  must  be  equal  to  the  work  done  on  the  piston  which  is 


where 


(A.3) 


But  from  the  shock  wave  equations  derived  from  the  condition  of  the 
Rankine-Bugoniot  equations  we  have 


and 


U 

u 


-  1 4  6  y 

'  5(jT-  1) 


( r  -  1.4) 


6  4  y _ 

5(y  -  l) 


(A.4) 
(A. 5) 


c  §(y  “  l) 

6  4  y 


(A.6) 
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Consequently  for  a  flat  top  shook  wave 


(A.  7) 


and  the  total  energy  of  a  shook  wave  of  arbitrary  shape  is  given  approxi¬ 
mately  within  an  error  of  2  percent,  providing  y  is  no  larger  than  3,  by 
the  following  equation 


PRESSURE  PH3TCT!IlfA  El  SHOCK  TUBE 
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